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1 Introduction 


We study the existence, stability and localization of soliton type solutions for 
the Non Linear Schrodinger Equation (briefly, NLSE) in the semiclassical limit 
(that is for h —>■ O'*"), for a bounded domain with Dirichlet boundary conditions. 
In our framework, the problem takes the form 

^ ^ C'i(R+,iLi(L!,C)) 

'il;{0,x) = (j}h{x), X e Q, ( 1 ) 


ij{t,x) = 0 on Rq X 

n C being open and bounded, N ^ 3, a > 0, where (j)h{x) G .ffo(fl) is 
a suitable initial datum, and V is an external potential. Conditions for the 
nonlinear term W and the potential V are to be precised and discussed in the 
following sections. 

The NLSE in the presence of a potential is largely present in literature. In 
particular, it has been extensively studied the effect of the potential V on the 
existence and the profile of a stationary solution, that is a solution of the form 
%j){t,x) = U{x)e~i‘^*', uj = A/2, where U solves the equation 

-h^AU + + 2V{x)U = XU. 

The first attempt to this direction is the work of Eloer and Weinstein [TS] , for 
the one dimensional cubic NLSE (with a generalization for higher dimensions 
and different nonlinearity in m) where, by means of a Lyapunov - Schmidt re¬ 
duction, it is proved that, if V has a non degenerate minimum, then a stationary 
solution exists, and this solution has a peak located at this minimum. Del, Pino 
and Felmer |14j showed that any (possibly degenerate) minimum of V generates 
a stationary solution. We also mention [3 HO], in which similar results are ob¬ 
tained with different techniques.Concerning global methods, in [24] Rabinowitz 
proved the existence of a stationary solution with a Mountain Pass argument. 
Later, Cingolani and Lazzo |13] proved that the Lusternik - Schnirelmann cate¬ 
gory of the minimal level of V gives a lower bound for the number of stationary 
solutions. The topological approach was also adopted in [5], where a more re¬ 
fined topological invariant is used, and in [^, where the presence of a negative 
potential allows the existence of a solution in the so called “zero mass” case. 

Another interesting feature is the influence of the domain in the stationary 
NLSE, when E = 0. In this case, a single-peaked solution can be constructed. 
In [^, Ni and Wei showed that the least energy solution for the equation 

— h'^Au + u = u^, u > 0 in D, (2) 

with homogeneous Dirichlet boundary condition, has a unique peak, located at 
a point Pfi with d{Ph,dil.) —>■ maxn d{P, 9D). Later, Wei [55] proved a result 
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that can be viewed as the converse of the forementioned theorem. Namely, the 
author showed that for any local maximum P of the distance from the boundary 
Sn, one can construct a single-peaked solution of (I2|) whose peak tends to P as 
/i —>■ 0+. The profile of the solution is, up to rescaling, close to the profile of the 
ground state solution of the limit problem 

—Au + u = u^, u > 0, 

in the whole of We also mention m in which the existence of a multi- 
peaked solution of dD) is proved. 

In the present work we follow a different approach, incorporating and ex¬ 
ploiting ideas found in II ism], where the problem (ED had been studied for 
the whole of for both cases: V = 0 (existence and stability), and V ^ 0 
(existence, stability and dynamics). 

We want to point out two main differences between this paper and [1]. First, 
dealing with a bounded domain and a bounded nonlinearity gives us enough 
compactness to easily prove an orbital stability result. In particular our result 
is also true for positive nonlinearities, which is forbidden in the whole space by 
Pohozaev’s Theorem. 

This difference is less evident in what comes after, since, when dealing with 
the semiclassical limit, we have to face a limit problem in the whole space K^, 
which is the same of [5], so we have to reintroduce some hypothesis. 

However, our orbital stability result could be used when studying other 
situations, e.g. fixing h and looking for solitons of prescribed, possibly large, 

norm (the so called large solitons). If one is interested in these topics, we 
recommend the nice paper by Noris, Tavares, Verzini [23], which deals with 
orbital stability for solitons with prescribed norms in bounded domains, and 
the references therein. We point out that in |23) . the authors work with pure 
power nonlinearities, so they can have some compactness loss -for L^-critical and 
supercritical powers- even if the domain is bounded. In particular, for critical 
and supercritical powers they have orbital stability only for small norms. 

The second point in which the bounded domain marks a difference with the 
paper [i], is pointed out in the Appendix. In fact, when trying to describe dy¬ 
namics, immediately it appears a repulsive effect of the boundary. The soliton 
is affected by a force oriented with the inward normal to the boundary. Unfor¬ 
tunately, the value of the force depends on the norm of the solution on the 
boundary, so, at the moment, we were not able to give a quantitative estimate 
of this repulsive force, and further efforts are needed. 

According to this line of thought, we have divided the present work into 
three sections and an appendix: 

In Section 2, existence and orbital stability results are obtained for the case 
V = 0, by referring to the related eigenvalue problem 


-AU + W\U) = XU, inn 


U = 0, on dn, 


(3) 
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given that a solution U{x) in i?o(n) of ([3]) results to a solution i/) = U{x)e~^ 
of (H]), with initial condition 'ip{0,x) = U(x). These results are summarized in 
Proposition 6. One should notice that the relative proofs work without having 
to impose the usual restriction 2<p<2+^, by relaxing the restriction to 
2 < p < 2* = instead. It is the boundedness of the domain that allows us 
in doing so. 

In Section 3, where we assume the presence of an external potential, our basic 
result, obtained by means of a rescaling procedure, is to prove localization in 
the sense that if we start with an initial datum close to a ground state solution 
U of 


- hAC/ + 


/i“+i 


W'{U) 


LO 

/l“+l 


U, in n 


U = 0, on dfl, 


(4) 


the corresponding solution of m will keep its profile along the motion, pro¬ 
vided that h is sufficiently small. Here and in what follows, the restriction 
2<p<2-|-^ is imposed, since we need to face a limit problem in 

In Section 4, an modular localization result is obtained for the case 
V ^ 0, and for both cases: the unbounded and the bounded one. When we 
work on the whole of we start with a ground state solution Ui of the 
counterpart of 0, proving that a solution of the counterpart of o, with 
initial condition close to C/i, preserves its basic modular profile as time 
passes, in the sense that given e > 0, for all t > 0, the ratio of the squared 
norm of \'S/uh{t,x)\ with respect to the complement of a suitable open ball 
over the squared norm of \Vufi{t,x)\ with respect to the whole of is less 
than e for h sufficiently small, where Uh{t,x) is taken by the polar expression 
of namely tp{t,x) = Uhit,x)e''^^^*’^'>. The bounded case is treated by 

exploiting ideas developed for the problem (Section 3). 

Finally, and as we mentioned before, in the Appedix it is described an at¬ 
tempt to study dynamics in the frame of a bounded domain, where we encoun¬ 
tered difficulties due to computational complications related to the action of 
W on the motion as well as to the repulsive effect of the boundary. 


2 The case V = 0 

2.1 Existence 

For simplicity of the exposition we assume h = 1. As it has been already said, 
the case H = 0 is related to problem (|3]), and a solution u(x) in ffgifi) of ([3|) 
results to a solution ^ = u{x)e~~ of 0 with initial condition ^/’(0,a:) = u{x). 
Notice that a minimizer of 

J{u)= J (^^\Vu\^ + W{u)^dx 
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on S'cr = |u S HQ{rt) : ||M||i 2 (Q) = cr|, for some fixed cr > 0, is a solution of ([3]), 
for suitable A. Thus we focus on the existence of such a minimizer. We impose 
on W the following conditions: 

Condition 1 W is a , bounded and even map ffi. — R. 

Condition 2 |lT'(s)| < c\sY’~^, 2 < p < 2* = where c is a suitable 

positive constant. 

Remark 3 Although it is intuitively quite clear the construction of such maps, 
an easy concrete example is furnished by choosing W = sin(s^), for p > 0, and 
evenly expanding it on the whole of R. We also stress the fact that a bounded 
W is related to the global well posedness results by Cazenave. 

Notice that 


— oo < pL = inf J(u). 


(5) 


If {un} is a minimizing sequence in So- for J{u), that is J(un) —>■ /r, it is evident 
that {un} is bounded in thus, up to a subsequence, Un ^ u € 

and —>■ u in The latter implies ||un|l 2 , 2 (n) ^ ll^llL2(a)’ ^ ^ 

Next, we obtain a similar result to Proposition 11 in [3], 

Proposition 4 //{?«„} is a minimizing sequence in So- for J, that is J{wn) —t 
pL, satisfying the constrained P - S condition, that is, there exists a real sequence 
Xn of Lagrange multipliers such that 


AWn + W'{Wn) - XnWn = Cr„ 0, 


( 6 ) 


then Xn is bounded. 


Proof. Since, as we saw above, Wn is bounded in dH]) implies 



where by H-H^ is denoted the dual norm for Hq{LI). We have 




\VWnf + 2W{Wn) - 2W{Wn) + W'{Wn)Wn - A„W^) dx 


2J(w„)-A„CT^+ / {W'{Wn)Wn - 2W{Wn)) dx 0. 


n 
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Notice that J{wn) is bounded, and because of Condition 2, 


{W'{Wn)Wn - ‘2W{Wn)) dx 


< 



< 


Cl ||wn||^i(Q) + 2fcmeas (fi) < +oo, 


where k is an upper bound of \W\. Thus A„ is bounded. ■ 

By Ekeland’s principle, if {ii„} is a minimizing sequence in for J(u), we 
may assume that it satisfies the constrained P - S condition, that is, there exists 
a real sequence A„ so that ([H]) holds. Because of Proposition 4, A„ is bounded, 
and the following hold 


Art A 

Un u in 

Un —>■ 11 in L^{VL) for 1 < p < 2*. 


We have already shown that u G Scr- Thus, u 0. Next, we show that 


Am + W'{u) = Xu. 


(7) 


To this end, if p is a test function, combining the three considerations above 
with Condition 2, we have 



implying ©. 


Notice next that due to Condition 2, the Nemytskii operator 

(fl), 2<t<2*, 

is continuous, whereas Un ^ u in L*(n), for 2 < t < 2*. Thus, 



proving that J{u) = /i. 

This completes the proof for the existence of a non trivial solution of (l3|), 
for suitable A. In fact, the weak convergence m„ ^ m turns out to be a strong 
one: Since J(m„) — J(m), f^(W(un) — W(u))dx —>■ 0, we obtain ||Mn||//i(Q) —>■ 
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||m||^i(q), thus proving that —>■ m in iJg(ri). Since W has been assumed even, 

we may take a non trivial nonnegative solution of (jS]) . By Harnack’s inequality, 
this solution is strictly positive on We thus obtain a positive solution u £ 
for problem ([3]), for suitable A. The wave function 

•4)(t,x) = u{x)e , w = A/2 (8) 

is a stationary solution of o, for h = 1, V = 0, with initial condition (j){x) = 
ijj{0,x) = u{x). Evidently, —w = A/2, is a stationary solution of (P), 
too. 

2.2 Stability 

We turn next our attention to the stability of the stationary solution. To this 
end, we focus on the reduced form of 

2^^ = in R+ X O, 

^/>(0,a;) = (/(a;), (9) 

x) = 0 on Mq X dfl, 

by taking, as it was mentioned above, h = 1. The different time slices iptix) 
of each solution of ([9]), where such a solution may be understood as the time 
evolution of some initial condition ipo{x), could be thought of as elements of a 
proper phase space X C L^(r2,C), with the set 

r = Irt(x)e*®, 9 £ R/27rZ, u £ Sa, J{u) = ^ = inf J{w) 1 (10) 

( weS„ j 

being an invariant (under evolution) manifold of X. Evidently, ±?l(x) £ T. 

To make the description of all this more clear, one should notice that if 
ipto (x) is a time slice of a solution x) of (P), the evolution map is defined by 

UtiJtoix) = tpto+t{x), 

meaning that this time slice might be considered as the initial condition of the 
solution tpi{t,x) = tp{t + to,x). Now, if u{x)e^^ £ T, then u is a solution of (P, 
with suitable A, and, at the same time, u(x)e*® is the initial condition of the 
solution '!/(t,x) = it(x)e*(^“'^‘/^) of (P. Since M(x)e*(®“^‘/^^ G T for each t ^ 0, 
the invariance of T follows. We are going to prove orbital stability of u(x)e , 
w = A/2, following the definition of orbital stability found in [12], meaning that 
r is stable in the following sense: 

Ve > 0, 3 (5 > 0 such that if ip{t, x) is a solution of P satisfying 


inf ||?/(0, 


lm(a) < then Vt A Oinf Wipit^x) - 


(n) 


< e. 


( 11 ) 
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The name orbital could be misleading in the present framework: our non¬ 
linearity does not ensure the uniqueness of the ground state. Moreover, the 
problem is not invariant under translations. Thus it could be that the set T 
reduces to a single function or to a discrete set of ground states (up to multi¬ 
plication by a unitary complex number). Anyhow, we follow the approach by 
fn\ and by [1], so we will keep the term orbital stability. In Lemmas 21 and 22, 
we will see how the orbital stability result leads to a precise description of the 
modulus of a solution tjj{t,x) starting from a suitable initial datum. 

Notice that T is bounded in since for each of its elements w = 

w{x)e^^, w{x) is a constrained minimizer of J, whereas W is bounded. Notice 
that we may take w{x) > 0. 

Suppose r is not stable. Then 3 e > 0, and sequences —)• 0+, ipnit^x) of 

solutions of dH), and > 0 such that 

inf ^„(0,a ;)ini Uriitn,x) - >£■ (12) 

Notice that the first inequality of m implies 

inf ||V'„(0,a:) - w|| < C'(5„ ^ 0, 

where C is the Sobolev constant satisfying ||•||^ 2 (f^) < Thus, we 

may obtain a sequence Wn in T, such that 

||?/’n(0,x)-w„||^^^^^ ^>0. (13) 

We express now ipn(t,x) in polar form, namely = Un{t, with 

Unit,x) = \'ip„(t,x)\, \/t A 0. Since ||wra|| = cr, (IT^ implies that u„(0,a:) is 

(H) 

bounded in and at least up to a subsequence, still denoted by u„(0,x), 

||wri(0, a:)|| —>• M > 0. Rewriting (|T^ in its squared form, and taking into 

L^(n) 

consideration that 


Un{0: x)\Wn{x)\dx < ||Mr!,(0,x) 


T2(n) 


we take 

0 > - 2Mcr + a^ = {M - af, 

thus obtaining M = a. 

The polar form = u{t, turns (|9]) into the system 


Am 


W'{u) 

2 


|+i|v,n.=o 


dtu^ + V • (u^Vs) = 0, in Kg x Q, 


M(0,x)e*®^°’“^ = (j){x), u{t,x) = 0 on Kg x dO., 


(14) 
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with the two equations of m being the Euler - Lagrange equations of the action 
functional 


A{u,s) = - jj \Vu\^ dxdt + ^ JJ W{u)dxdt + ^ jj 

The total energy is given by 


ds 




u^dxdt. 

(15) 



E{i^) = E{u, s) = ^ |Vu|" + 1 Vs|' + W{u)^ dx, 

(16) 

that is. 

E{iIj) = E{u,s) = J{u) + - f u^\Vsfdx 

2 Jn 

(17) 

Independence of time for the energy and for the charge imply that for a solution 
= u{t,x)E'‘^^*''^'> of dni), it holds 


d f 

— / u{t, x)‘^dx = 0 
dt Jq 

(18) 


|i5(„,s) = 0. 

(19) 

Equivalently, (ITOl) and (ITOl) can be expressed as 




(20) 


E{'ip{t,x)) = E{(j){x)) 

(21) 


for all t > 0. Noteworthy, for stationary solution, m yields E{tp) = J{u). 

Returning to the sequence ipnit^x) satisfying (fT^ . we may assume, as we 
saw, that ||u„(0,a;)|| —>■ a, that is ||it„(t,a;)|| ^ a, for f > 0, because of 

(f2) L^{d) 

([201). We want to show that {««}„ is a minimizing sequence for the functional 
J on the constraint ||u|| = cr. 

One should notice that the first inequality of m, combined with the bound¬ 
edness of r ensure that ipn{0,x) is bounded in TLq)!!). Since W is bounded, 
CZl) ensures that i?('0„(O,a;)) is bounded, and because of (1^ . E{'4>n{t,x)) is 
bounded, for all n and all f > 0. In particular, E{ipn{tn,x)) is bounded. A 
new application of (II3, ensures now that Un{tn,x) is bounded in E[q{^). The 
sequence Un(tn,x) = where = n., *1,11 -, is in Sa. We 


i (^n 1^) 


L^{n) 


have, writing for simplicity Unj Un instead of Un{tn^^)i respectively, 
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for suitable In = ln{x) S (0, 1), and because of Condition 2, 

I J(u„) - J(u„)| < ^\an-l\ [ \yunf dx+ [ \W{Un)-W{Un)\dx 

^ JQ Jo, 

= \Wl-M [ \^Un\^ dx 

J Q 

\^n (^n^n ^n)^n)| 

Jn 

+ \o^n ~ ^\!"n (X ~ ^n) Q^n] 


^-^\Un\^dx 


0 , 


since in the right hand side of the last inequality, the two summands are products 
of a zero sequence by a bounded one. Thus, J{un) — J(un) 0. We return 
now to 

||V'„(0,cr)-rh„|| ^0, (22) 

Hq (Si) 

which, as a result of the triangle inequality combined with the boundedness of 
||V'n(0,a;)|| +||w„(x)|| , readily gives 

Hl(Si) Hl(Si) 

||'0n(O,a:)||^ - ||w„(a;)|l^ 0, 


h1(S2) 


h1{S2) 


that is, 


|Vu„(0,a;)| +it^(0,a;)|Vs„(0,a;)| -|Vw„(x)| da; 0. (23) 


Jo 

We claim that 


M^(0, x) |Vs„(0, x)\^ dx —>■ 0. 


If not so, up to a subsequence. 


lo 


|Vu„(0,a;)|"-|Vu>„(a;)r 


dx ^ k < 0. 


(24) 


(25) 


Combining convergence of Un(0,x) — Wn(x) to 0, with Condition 2, we have 


/ [W(u„(0,a;)) - IT (?«„(a;))] 0. 

Jo 

Now (1^ and ((^ give 

J (Un(0, x)) — J {Wn{x)) —)> k/2 < 0. 


(26) 


(27) 
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However, as we have shown above, 


thus obtaining 


J(w„(0,a;)) - J{un{0,x)) 0, 


J{Un{0,x)) - J {Wn{x)) k/2 < 0, 


(28) 

(29) 


an absurdity, since Mn(0,x) S Sa, and Wn{x) is a Sa minimizer of J. Thus 
holds, and because of (1^ and (1^ . we get 


J('U„(0,x)) - 5 > fi. 

£’(V'„(0, x)) - E{wn) 0, 

x)) - E{Wn) 0. 


(30) 


Thus, 

and by m we have 
From (1501) we obtain 

/ u^(t„,x) |Vs„(tn,x)|^dx0. (31) 

Jn 

To see this, let us assume that this is not the case. Then, up to a subsequence, 

/ u^(t„,x) |Vs„(<„,x)|^(ixp > 0. 

Jn 

Then (l30)) implies that 

J (■u„(<„,x)) - J (w„(x)) -p < 0. 

Since 


'^(^n(^ni^)) Jij^nijn^xy) y 0, 


we obtain 


x')') J (^Wn{x)) y p <C 0, 

which is absurd, for the same reason as above. That is, (IQTl) holds, resulting to 

J (Un') — J(Un{tniX^^ p. 


thus implying 


Ji^nijn ^ x')') y p. 


In other words, we may consider u„ = Un{tn,x) as being a minimizing se¬ 
quence in Sa- for J{u). As such, by the previous discussion, Un(tn,x) —>■ 
u'(x), with u'(x) € Sa, being a minimizer of J{u). Now (I31|) ensures that 
||'0n(in,a^)ll^fi(n) ll'“^(^)lliri(n)! by lower semicontinuity of the norm we 
finally obtain that il)n{tn,x) —>■ u'{x), thus proving orbital stability of the sta¬ 
tionary solution. 
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Remark 5 As we claimed before, once '0(0, a;) is sufficiently close to F, ip(t,x), 
i ^ 0, are all close to F. This implies that the modulus \ip(t, x)\ is, for any t ^ 0, 
close to a (possibly different) ground state of the NLSE. This consequence is the 
key tool for proving two main results of the next section. Lemmas 21 and 22. 
Thus, the orbital stability in this context becomes a localization result for the 
moduli \'if{t,x)\. 


We summarize the existence and stability results in the following: 


Proposition 6 The problem 

= -Y** + ■" < X n, 


iIj{0,x) = (j)h{x), X G LI, (32) 

tp{t,x) = 0 on Rq X dLl, 

where h > 0, and (fhix) being a suitable initial datum, admits a stationary 
solution flit, x) of the form Uh{x)e~^^*. More concretely, Uh{x) is obtained as a 
solution of the eigenvalue problem 


h ^ 1 

2^“^ 2/i“+i 


W'(u) = ku. 


LI 


u = 0, on dLl, 


(33) 


for suitable k. In addition, f(t,x) is stable in the sense of RJ|) . 


We give next the definition of a solitary wave with respect to a bounded 
domain Ll. To this end, we have to define first the notion of the barycenter 
of a family of states ft{x), t ^ 0, whose members are obtained by the ’’time” 
evolution of an initial state foix), in the frame of a proper phase space X C 
L‘^{Ll,<C). 

Definition 7 For ft(x), t ^ 0, as above, its barycenter, q{t), is defined by the 
relation 

,(() = (34) 

In\Mx)\ dx 

Remark 8 An analogous definition of the barycenter is given in under the 
condition that it makes sense. In our case, the definition of q(t) makes always 
sense, because Ll is bounded. 


Remark 9 Notice that q(t) does not belong to Ll necessarily, unless Ll has spe¬ 
cific geometric features. For instance, convexity of Ll would ensure that q(t) G Ll 
for all t ^ 0. 
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Definition 10 The state tp = ipoix) in the phase space X C L‘^(fl,C), is called 
a ’’solitary wave” in the frame of a dynamical system Utip = ''Ptix), t ^ 0, 
where U : x X ^ X is the evolution map, if: Given e > 0, we may find 

k{e) > 0 such that for each t ^ 0, there exists a neighborhood Ve,t of q(t) with 
meas[Tl — fl O)] ^ k{e), and 

[ \ipt{x)\^ dx - f \pJtix)f dx < e. (35) 

JQ Jve,tr\n 

Remark 11 It is easy to see that the stationary solution of iStH) is a solitary 
wave in the above sense: The barycenter in this case is fixed for all t, and 
one needs to suitably blow up a given neighborhood of it, in order to meet the 
requirements of the above definition. 

3 localization for V 0 


For the rest of the exposition, we assume without loss of generality that 0 € fl. 
We also restrict p in Condition 2 so that 2 < p < 2 + , and we impose on W 

the additional condition: 

Condition 12 3so such that lF(so ) < 0, 

and on V the following one: 

Condition 13 V{x) € C^{fl), is nonnegative. 

Remark 14 With these restrictions and additional conditions on the nonlin¬ 
earity and the potential, problem {Ip is globally well posed (see \11\. Thm 3.3.1, 
and Thm 3.4.1}), with the energy and the mass remaining constant in time. 
Commenting especially on Condition 12, we notice that in the opposite case, 
there would exist no nontrvial stationary solution for the related K.^ problem, 
undermining thus the basic tool we will use in order to prove the main stability 
result. In the previous section, where we had posed weaker conditions on fl, 
W, and p, we were led to nontrivial stationary solution without any particular 
comment on X. It is trivial to see that for W having positive lower bound, A is 
positive, too. Under the additional conditions imposed above on fl, W, and p, 
we can prove, as we will see soon, that we may obtain a nontrvial solution of 
0) sitting in So- for suitable a, with X = X{a) < 0, provided that fl contains a 
suitably big open ball B{0,r{a)) centered at 0 with radius r{a). 

3.1 Rescalings 

We set /3 = 1 + ^. For h < 1, we define the inflated domain 

flh = {xGR^ : h^x e fl} . 
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If u is a solution of the stationary problem 


— Am + II^'(m) = WM, in 

M = 0, on drih, 

then Vh{x) = v{^) is a solution of the stationary problem 

— hAu + W'{u) = ^ M, in 11 

ha+l ^ > /ja+1 

u = 0, on dVl. 


(36) 


(37) 


Furthermore, we define the functionals 

^h{u) = J u'^{x)dx, Cn^{u) = J v?{x)dx, 

dx, Jq^{u) = / 
Jq 


Jh{u) = 




^2 

y |Vm|" + IF4m) 


^\Vut + W{u) 


dx. 


where Wh{u) = -^W{u). 

We have the following identities: 

Jh(vh) = h~°‘Jn^{v), Ch{vh) = Cvihix). 

We next define 

m{h, f2) := inf Jh- 
Ch—'^ 


(38) 


Lemma 15 Forfli,fl 2 two bounded domains as described above, withfli C 122 , 
it holds m{h, 122 ) A rn{h, 12i). 

Proof. It is straightforward, since i7Q(12i) C i7Q(122 ) by extending a function 
in i7Q(12i) into a function in Hq(P ,2 ) by zero on 122 \12i- ■ 

Lemma 16 m[h,Ft) = h~°‘m{l,flh), where to(1,12;i) = inf 
Proof. By rescaling. ■ 

Notice that the conditions satisfied by p, W correspond to the prerequisites 
for the existence result given in [4] to hold. Namely: 

Lemma 17 There exists some a > 0 such that for all a > a a positive mini- 
mizer Ua G exists for J{u) over all u G ||m||^ 2 (rw) = cr. In 

faet, 11 Wo-11 2,2 (RiV) = <J, J{ua) < 0, and Ua is a solution of the version of (0), 
with A < 0. 
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We consider an increasing sequence {r„}, with r„ —>■ cxd. Let 

mcr{B{0,rn)) ■■= inf J{w). 

ll-llr2(B(0.r„))=- 

Trivially, as Lemma 15 indicates, we see that 

m^{B{0,ri)) ^ ma{B{0,r2)) ^ ^ rriai^^) = J{w) >—oo. (39) 

ll«'lli2(RiV)=<T 

Lemma 18 lim m^(B{0^rn)) = ma-(JS.^). 

n—¥oo 

Proof. By Lemma 17, mcr(R^) is attained by some u G 

with ||w|| 2 , 2 (rW) = cr. Actually u is radial ([H], Thm. II.l and Rem. 11.3). 
For each n, we may choose a C°°(M^) real function Xn satisfying 

Xn = 1 if |a;| < rn/2, 

Xn = 0 if |a;| ^ r„, (40) 

|Vx„| < 4/r„. 


We define Wn = XnU- For suitable tn > 0, we have ||in'U'„||j;^ 2 (]aiV) = 
||inWri||i2(B(0,r„)) = Setting Un = tnWn, dM]) yields 

J{un) ^ m^{B{0,rn)) ^ m^(K^). (41) 

We want to prove lim J{un) = that will finish the proof. We have 

n—¥oo 

Wn ^ u in L^(K^). Thus —>■ 1, and Un ^ u in L^(R^). We also have 


/ |Vw„-Vu|^(ix = / |(Vx„) u + (xn - 1 )Vm|^ dx 

< 2 / |(Vxn)M|^ dx + 2 / |(x„ - l)v?l|^ dx 

JR'^ jR" 

< 32/r^ f |?l|^dx + 2 f |Vu|^dx—>-0 

J\x\>rn/2 J\x\>rn/2 

as n ^ 00 . So Vwn Vu in L^, and Vu„ —?> Vu in too. Thus tt„ —^ ?I in 

This is combined with the continuity of the Nemytskii operator 


to ensure that 


W : L* (R^) (R^) , 2 < t < 2*, 


/ VF(u„)dx—>■ / kF(u)dx. 

Im.^ jR'v 


Thus lim J{un) = m„{R^), and the proof has been completed. 


Remark 19 The above Lemma makes clear the final assertion of Rem. 1^: If 
Q contains a suitably big open ball B(0, r(a)) with W < a, then TOo-( 1, fl) has to 
be negative, and so has to be the eigenvalue A = A(cr) related to 
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Lemma 20 For a sequence of positive numbers hk —^ 0, for k —> oo, it holds 

lim 

k—^oo 

where ini Jq^ 

Proof. Combine Lemmas 15 and 18. ■ 


3.2 localization 

To facilitate exposition, we make the harmless assumption that a = 1, thus 
suppressing subindices in all involved infima m. We have the following: 

Lemma 21 For any e > 0, there exist 6 = 6{e), ho = /io(e) > 0, and R = 
R{e) > 0 such that, for any 0 < h < ho{e), there is an open ball B{qh,hl^R) C 
n so that for any u S i7g(n) with Ch{u) = 1, and Jh{u) < m{h,fl) + 6h~°‘, 

~UN~8 I ^ ^ ^ 

h Jn\B{q^,hPR) 

to hold. 


Proof. We argue by contradiction. Assuming the contrary, there exists £ > 0 
such that for any r > 0 we may find sequences (5„ = 6n{r), hn = hn{r) —>• O’*', and 
a sequence Uh„ = Uh„{r) € Hq{D,) with Ch„iuh„) = 1, Jh„iuh„) < + 

Snh~°' such that, for all open balls B{qr, h^r) C 11, 

Uh„dx > sh^d 

r.hnr) 

to hold. For each n, we now pass to the counterpart of Uh,^, denoted 
by Un, that is, Un{x) = Uh„{h^x). Combining (j38p and Lemma 16, we have 
Cnh„ (un) = 1, JQhn (Wn) < "i(l> and 


ln\B{q 


/ ul^dx > e, (42) 

According to Lemma 15 in [3, there is a d > 0, and an open ball B{q,R') in 
such that, for each w in with ||cc||i 2 (RiV) = 1, Jiw) < m(M^) + (5, 


/ w^dx < e 

Js.^\B{q,R') 

to hold. Take now d„(i?'), hn{R') as above. Because of Lemma 20, for n big 
enough, we may ensure that B(q,R') C and < 

to(R^) + S. Extending by 0 outside i}h„(R'), we obtain a F[^(R^) function 
meeting the requirements of Lemma 15 in [7], thus 

/ u^dx = / u^dx < e, 

JR^\Biq,R') 
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that contradicts (H^ . ■ 

Lemma 21 makes obvious the following: 


Lemma 22 For any e > 0, there exist S = (5(e), hg = hQ{e) > 0, and R = 
R{e) > 0 such that, for any 0 < h < ho{e), there is an open ball B{qh,hl^R) C 
0 so that for a solution ip(t,x) of (Q]) with Ch{\'4’(t,x)\) = 1, and Jh{\'f’{t,x)\) < 
m{h, n) + 6h~°^, for each t G RJ, 


1 


lQ.\B(qh,hf>R) 


|'0(t, x)\^ dx < £ 


to hold. 


The correlation of the solutions of the equations (l36l) and (l37)) . combined 
with (1551) . ensure the existence of ”ground state” solutions of (EZl), that is, 
solutions that are minimizers of Jh{,u), with u satisfying Ch{u) = 1. 

We define next, the following set of admissible initial data, for given K, 
h>G: 

with Ufi(0,x) = (U + w)(x) 

t/ is a ground state solution of (1^ . and w G H^(n) s.t. 

||Vw(0,i)||i. <Kh"fl^ 

' {x)u\{0,x)dx < Kh^^ 

We prove next the basic stability result. 

Proposition 23 Given e > 0, there exists ho = / 10 (e) > 0, and R = i?(e) > 0 
such that, for any 0 < h < ho{£), there is an open ball Bfcfh, h^R) C ft so 
that for a solution 'fi{t,x) of (QP with Ch{\ip(t,x)\) = 1, and with initial data 
ijj{0,x) G Bf^, where K is a positive fixed number, it holds 

-rm f \f;{t,x)\^dx <£, 

h P jQ.\B(qu,hPR) 

for any t G . 

Proof. Because of conservation of energy, we have 


E{ip{t,x)) = 


< 


< 


E{tlj{0,x)) 

h^^Jh{uh{Q,x)) 



|Vsh(0,a:)|" 

2 


+ V{x) 


dx 


h^f^Jh{uh(f),x)) + ^h^P + Kh^P 

h^^Jh{uh{Q,x))+Ch^^ 
h^l^Jh{U + w)+Ch^^ 
h^^{m{h,Vt)+C'Kh°‘+ C) 
h^f^ [mih, n) + h-°‘ {h'^°‘C'K + Ch°‘)] , 
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since a Mean Value Theorem application ensures that Jh{U + w) < m{h, Q) + 
C'Kh°‘, for suitably small h. More precisely, since Jh is , we may find some 
77 G (0,1) so that 


Jh{U+ w) - Jh{U) = J^(t7 + ryw) [w] 

= / V{U + 'qw)'S/wdx 


W'{U + r]w)wdx. 


We have 

V(C/ + r]w)Vwdx 


< livc/ll L2(n) 

(44) 

since < 1. Because of Condition 2, 


W'{U + r]w)wdx 


< C 2 ||w||^p(f2) < Ca . 


(45) 


Combining (l44l) and (l4^ . we obtain the desired inequality. Thus 


Jh{uh{t,x)) = h \E{il}{t,x)) - j ul{t,x) 


\Vsh{t,x)f 


■V(x) 


dx 


< [m{h,n)+h-°‘{h^°‘C'K+ Ch°‘)] , 


since V{x) > 0. We use now Lemma 22, by choosing Hq small enough in order 
to ensure K^°^C'K + Ch^ < (5(e), and the result follows. ■ 

Corollary 24 Given e > 0, if q{t) is the barycenter of a solution ifltyx) of (Qp 
with Ch{\il>{t, x)|) = 1, and with initial data '0(0, x) G , where K is a positive 
fixed number, then for any t G K))", q{t) G B(cfh, de + h^R), where h and R are 
as in Proposition 23, and d = diam{VL). 

Proof. One has for any t G M))", 

Tni, — rr 

\qh-q{t)\ = 


In {qh - x) \ip{t,x)f dx 

1 

II 

[ {qh - x)\ip{t,x)f dx 

fa \fj{t,x)f dx 


Jn 


< 


h [ \qh - x\\'ip{t,x)f dx 
Jn 

h-Np 


!n\B(qH,hf>H)\'^h-x\\'iP{t,x)\ dx 

+ L{qu,hf<R) \dh-x\ \'ll^{t,x)\^ dx 
< h-^^ [dh^^e + h>^Rh^^] =de + h^R, 


thus proving the claim. 
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Remark 25 The above Corollary may he interpreted as follows: If we choose 
a pretty small e > 0, and we make a consequent choice of small h, then 
B(qh,h^R) C Ah = B{qh,de + h^R) fl Q., with meas{fl — Ah) > 0. In other 
words, we have found a neighborhood of q(f), for any t S Kj, as in Def. 10, 
in the sense that a solution tfitjx) of ([I]j with a perturbed initial state, meeting 
the requirements of Proposition 23, concentrates basically on this neighborhood 
of q{t), exhibiting a behavior quite similar to a solitary wave. 


4 localization for V 0 

4.1 The case 12 = 

We start with the following assumption: The problem 

"'■W = -T^*+ • 

'ip{0,x) = (fhix), (46) 

\tf{t,x)\ —5> 0, as \x\ —>■ oo, 

considered as the analogue of © admits a unique solution 

gC°(R, H^{R^))nC\R,L^{R^)) (47) 

(see m, [TT] or [18] for sufficient conditions.) We also impose on W, V further 
conditions, namely: 

Condition 26 W is , with |W"| < c for some c> 0, 2 < p < 2 + -^ 
Condition 27 V : R^ R is a nonnegative function 


In order to proceed, we will need the following lemma: 


Lemma 28 For every e > 0, there exists R = R{e) > 0 such that for every 
ground state U, there exists q{U) G such that 




(|VC/|Vt/^) dx<e 


Proof. If we assume the contrary, then we may find £ > 0 such that we may 
have a sequence of pairs (i?„ > 0, t/„ ground state) so that for each q £ R^ 


lR^\B{q{U),R„) 


(jyUnf + U^'j dx>e, 
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thus obtaining 


inf [ 

jRN\B{q(U),R„) 


(|VC/„|" + dx > e. 


(48) 


Then {[/«}„ is a minimizing sequence, and by concentration compactness we 
know that {C/n}„ is relatively compact up to a translation by {<?«}„ G Thus 
there exists a ground state U with Un[ — 9n) —>■ C/ in and 


inf / 

jR’^\B{q(U),R^) 


V17„r + U^. dx 


< 


vunr + u^)dx 

{\yUnf + Ul){x - qn)dx 

[ (\VU\^ + uAdx + 0 r,{l)= 0 ^{l), 

Jr^\b{o,r„) ^ ' 


contradicting 


Lemma 29 For every £ > 0, there exist R = i?(e) > 0, <5 = (5(£) > 0 such that, 
for any u £ H S'o-, we can find a point q — q{u) € such that 


jR^\B{q,R) ^ ^ 

Irn |Vu|^ {x)dx 

where m = m(K^) (see Lemma 21), J™-+^ = g i7^(R^)/J(u) < m + S}. 

Proof. Exploiting Rel. (50) in [7], we obtain a point q = q{u) G K'^ and a 
radial ground state solution U such that 


u{x) = U{x — q) + w, with < C'£, (50) 

where C is a constant not depending on U. According to the previous Lemma, 
we may find R > 0 and a point q = q{U) such that 

[ (\VUf + U^)dx (51) 

jR^\B{q,R) ^ ^ Cj 

where ci is the Sobolev constant related to the embedding iJ^(K^) ^ 

If we choose R big enough, then B{q,R) C B{0,R), resulting to 


/r'v\B(o,R) 


dx 


< 


We have 

/rn\b(oA) (^^dx <^1 /rn\b(o,b) +U^dx 


/Riv |Vt/| {x)dx 


< 


U^dx 


< Cs. 


(52) 


(53) 
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Now 


i 


R'^\B(q,R) 


|Vm|^ {x)dx < 


R^\B(q,R) 


R^\B{q,R) 


\VUf {x - q)dx 



R^\B(0,R) 


\VU\^{x)dx 



By (l50l) . (l5^ and (l54l) . we get the claim. One should notice that R does not 
depend on u, U. ■ 

Our main objective in this subsection is to prove an modular stability 
result of the solution of (I46p with suitable initial data; more precisely, we prove 
that, for fixed t G Mj, this solution is a function on with one peak localized 
in a ball with its center depending on t whereas its radius not. To this end, it is 
sufficient to assume that ((46l) admits global solutions 'tp{t,x) G 
satisfying the conservation of the energy and of the norm. 

Given K > 0, h > 0, we define the following set of admissible data: 


with M,i(0, x) = (U + w) (^) 

q G 17 is a ground state solution, and w G i7^(R^) s.t. 
||t7 + w \\^2 = \\U \\^2 = cr, and HicH^i < 
i|Vs/j(0,x)||^o„ < K for all h 
/jjN V{x)u\{Q,x)dx < 


(55) 


We next study the rescaling properties of the internal energy 



and of the norm of a function u{x) having the form 



with = 1 + f. We have 
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and 


Jh{u) = 


rh2 


7r" 


^T.V 




dx 


^( 7 V- 2 )/ 3+2 


= h^P-' 


\V^v{0\^ + h^P-^W{v{0) 




7r" 


2 L 


|Vjz;(0|Vw^(^(0)l =h^^-^J{v). 


We can now describe the concentration properties of the modulus of the 
solution of (PSI) . 

Lemma 30 For any e > 0, there exist positive numbers 6 = (5(e), R = Ris) 
such that: for anyip(t,x) that solves with \ilj{t,h^x)\ € J"*+'*nS'cr, for all 

t, there exists a map qu '■ —>■ R'^ for which 


fgN \Vuh{t,x)\^ dx 


< e. 


(56) 


Proof. For fixed h and t, we set v{f) = \tp{t,h^f)\. By Lemma 29, there exist 
(5 > 0, i? > 0 and q = q{v) such that: if |^(f, h^x)| G j"^+^ n S^, then 


e > 




UlViv(Ol^df 

By a change of variables, we obtain 

/r"\B(9,_R) 1^5^ (^)I _ lR:^\B{hPq,hPR) ^)l ^2: 


e > 


> 


/ r(v |Vjn(0|^dC |Vu,,(t,x)|^dx 

/R'V\B(/tf'q^fl) ^)l dx 


/rn |Vu/,(t,x)| dx 


since h < 1. 


Setting qh{t) = h^q, we complete the proof. Notice that 
if, while R depends only on e. ■ 

Proposition 31 Let V G For every e > 0, there exists R > 0 and ho > 0 
such that, for any if{t,x) that solves with initial data if(0,x) G Brr, 

where h < ho, and for any t, there exists cfhit) G R^, for which 


\\Vuh{t,x )\\^2 JR^\B(qh{t)MR) 


\Vuh{t, x)\ dx < £. 


(57) 
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Proof. By the conservation law, the energy ^)) is constant with respect 

to t. Then we have 




< 


Eh{tjj{0,x)) 

Jh{uh{0,x)) 


/ ul(0,x) 

0 +V(x) 

/k" 

Z 


Muh{0,x)) + ^a^h^^ + Kh^^ 
hNP-°‘J{U + w) + h^>^C, 


where (7 is a suitable constant. By rescaling , and using that '0(0, x) G BK,h, and 
that ||w||^i < Kh°‘ implies J{U + w) < m + Kh°' (see the proof of Proposition 
23), we obtain 

Eh{ip{t,x)) = + w) + h^^C 

< (to + Kh°‘) + 

= (to + Kh°‘ + h°‘C) = (to + h°‘Ci ), 


where we have set Ci = K + C. Thus 
Jh{uh{t,x)) = Eh{tjj{t,x))- 


[ ul{t,x) 


/RN 

2 


dx 


since P > 0. By rescaling the inequality ([5^ . we get 

J {uh{t, h^x)) < TO + h°‘Ci. 

So, for h sufficiently small, we may apply Lemma 30, and get the claim. 


(58) 


4.2 The case f2 is bounded 

The case where is bounded is easily treated, once one makes the crucial remark 
that Lemma 29 has to replace Lemma 15 in [7], that it was used in the proof of 
Lemma 21. The rest of the proofs in the consequent Lemmas 32 and 33 and of 
the final Proposition 34 follow precisely the pattern of the proofs for Lemmas 
21, 22, and of Proposition 23, respectively. For completeness, we give below 
the precise statements, where we have assumed for simplicity, as in the case, 
that (7 = 1. 

Lemma 32 For any e > 0, there exist 6 = 6(e), hg = ho(e) > 0, and R = 
R{e) > 0 such that, for any 0 < h < ho(e), there is an open ball B(qh,hl^R) C 
so that for any u G Lfo(^) Ch(u) = 1, and Jh(u) < m(h,Vt) + 5h~°‘, 

In\B{qh,hPR) 

/nlVwl^da: 

to hold. 
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Lemma 33 For any e > 0, there exist S = (5(e), ho = /io(e) > 0, and R = 
R{e) > 0 such that, for any 0 < h < ho{e), there is an open ball B(qh,h^R) C 
fl so that for a solution ip{t,x) of (Q]) with Ch{\'4’it,x)\) = 1, and Jh{\'f’it,x)\) < 
m{h, 17) + 6h~°^, for each t S Kd, 

2 ;)| dx ^ ^ 

\ Vu{t, x)'^ dx 

to hold, where u(t,x) = \il){t,x)\. 

Proposition 34 Given e > 0, there exists ho = ho{e) > 0, and R = R{e) > 0 
such that, for any 0 < h < ho{e), there is an open ball B{(fh,h^R) C fl so 
that for a solution ifftjX) of (QP with C'ft,(|(/5(t, x)|) = 1, and with initial data 
ip{0,x) G , where K is a positive fixed number, it holds 

fn\B(qhMR) ^ ^ 

J^\\/u{t,x)f dx 

for any t G M()', where u{t,x) = \'tjj{t,x)\. 

5 Appendix 


In order to gain some control over the dynamics of the problem, that is, to 
try to formulate Newton’s equation describing the motion of the barycenter for 
a hxed value of h, one needs to express suitably q{t). To this end, a further 
assumption on W is made, namely that W{Q) = 0. For the sake of simplicity, 
we fix h = 1, a = 1 in what follows, the general case being straightforward. 
Also, we suppose that a solution tfitjX) is sufficiently smooth in order to make 
the integration by parts meaningful. Given this, the general case can be proved 
with minor technical efforts. Finally, we use the Einstein convention on the 
summation indices. 

We will use the Lagrangian formalism. Equation o is the Euler-Lagrange 
equation relative to the following Lagrangian density C: 

C = Reiiffdtfj) - i \Vfjf - W{\fj\) - V{x) 

By Noether’s theorem, there are continuity equations related to £, which we 
will use to derive an equation for the motion. In particular, we are interested 
in the following continuity equations: 


^|V’(t,x)p = -V-Im (V'VV’)) 

(59) 

^Im (^VV()) = -IV’PVE - V • T, 

(60) 
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where T is the so called energy stress tensor and has the form 


Re + - |VV>|" - -W' (1^1) \i,\ + WM) 


'^jk — R-S Sjk 

For an introduction to the Lagrangian formalism for equation m and continuity 
equations we refer to OlHllle]. 

In the light of equation [SHI and by divergence theorem, one has for j = 

= 4 / Xj\tp{t,x)f dx = [ Xjdt{\ip\'^)dx 
t/ o r2 

= — XjV ■ Im {t/jV-tp) dx 
JQ. 

= -[ V ■ [xjim (V'VV’)] + [ lm{'ipdx^'ip) 

Jn Jn 

= / Im {pjdxjtp) , 

Jn 

since ip(t,x) = 0 on dVl. Thus we have the momentum law 


q{t) = Im ( / tpVipdx^ . 


(61) 


For the second derivative of the center of mass, we have, by dsni) and by 
divergence theorem, 

qAt) =4- [ lm{AdxjAdx) = - [ dxi,Tjk{t,x)dx - [ x)\'^dxjV{x)dx 
™ Jn Jn Jn 

= [ Tjk{t,x) ■ Ukda - [ \ip{t,x)\'^dxjV{x)dx := h + h, 

Jon Jn 

n being the inward normal to 911. 

Let us use the polar form ipAyx) = u(t,x)e^^^*’^\ Then 

l2 = - [ u^dxj V(x)dx. 

Jn 


This appears to be a force term depending on the potential V. This, when the 
concentration parameter h —>■ 0, gives us the Newtonian law for the motion of 
a particle (see [7], where this approach is used in the whole space R^). 

Since u = 0 on the boundary (and since IT(0) = 0), the expression of T is 
simplified and the term Ii becomes 


h 



Tjk'^kdo' 





SjkA (u^) 


Ukda 




Again, because u = 0 on the boundary, by implicit function theorem, we have 
that Vu is orthogonal to d^l. In addition, we have by deflntion it = |'0| > 0, 
so whenever Vu 7^ 0, the inward pointing normal vector can be written as 


n = Thus 
|v«| 


= I {dxjU\Vu\-l-dx.u\Vu\)da =\- [ d:^^u\Vu\da = 
Jon ^ ^ Jdn 


'an 


iVitj^n^ticr. 


Concluding, we have 


q[t) = — f u'^S/V{x)dx + ^ f |Vitpndcr. (62) 

JQ 2 Jqq 

In the case of a bounded domain with Dirichlet boundary condition, it appears 
an extra term, which represents the centripetal force. Unfortunately, there are 
some obvious computational challenges concerning the last integral of (l62l) , and 
we cannot give a simple expression of this term, when h —>■ 0. As it was said in 
the Introduction, these challenges call for further work on the dynamics of the 
solution of O- 
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